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Abstract 

Recent anomalies in cosmic rays data, namely, from the PAMELA Collaboration, can be interpreted in terms 
of TeV scale decaying/annihilating dark matter. We analyze the impact of radiative corrections coming from the 
electroweak sector of the standard model on the spectrum of the final products at the interaction point. As an 
example, we consider virtual one loop corrections and real gauge bosons emission in the case of a very heavy vector 
boson annihilating into fermions. We find electroweak corrections that are relevant, but not as big as sometimes 
found in the literature; we relate this mismatch to the issue of gauge invariance. At scales much higher than the 
symmetry breaking scale, one loop electroweak effects are so big that eventually higher orders/resummations have 
to be considered: we advocate for the inclusion of these effects in parton shower Monte Carlo models aiming at the 
description of TeV scale physics. 

1 Introduction 

At the TeV scale and beyond, electroweak (EW) radiative corrections enter into the realm of nonperturbativity: one 
loop corrections relevant for the LHC can reach the 40 % level (see for instance Q]). It is surprising that the same 
electroweak radiative corrections produce small effects (typically less than 1 %) at LEP that probe the characteristic 
scale of the theory of 100 GeV and become huge at energies only 1-order of magnitude bigger. The reason for this is 
the presence of energy-growing contributions that, as has been pinpointed in [2J, are related to the infrared structure 
of the theory. More precisely, one loop corrections feature double logarithmic contributions oc aw log (i/s/mw), 
with yfs being the typical cm. energy of the process considered, and mw the weak scale of the order of W and Z 
gauge bosons masses; the weak scale itself acts in this case as an infrared regulator. Various interesting features of 
electroweak radiative corrections at energies much higher than the weak scale have been studied in the last ten years: 
noncancellation between real and virtual contributions, which is a unique feature of weak interactions [3], resummation 
of leading effects [4], relevance for phenomenology, and, in particular, for LHC processes [5]. 

Recent excesses observed by PAMELA [6], FERMI [7], and ATIC [8] (see also [9]) can be interpreted in terms 
of heavy- mass (1 TeV or more) dark matter (DM) annihilation or decay [TO]. Clearly, even if the final products are 
initially constituted by, say, an electron/positron pair sharing half of the cm. energy each, radiative virtual corrections 
and emission of additional particles in the final state will alter the injection spectrum at the interaction point. Then, 
the following relevant question arises. 

Assuming that physics below the DM mass is the standard model (SM) one, and assuming that the primary 
annihilation/decay process is known, what is the final products spectrum? 

Even if the physics describing the process is assumed to be perfectly known, the answer to this question is by 
no means trivial. The usual approach in the literature is to describe the effect of QCD and QED through analytical 
calculations and Monte Carlo generators like PYTHIA [11]; radiative corrections due to weak gauge bosons are usually 
neglected. However, including electroweak effects is important for at least two reasons. Qualitatively, since all SM particles are 
charged under the SU(2)l <8>U(l)y group, because of particle radiation the final spectrum will be composed of all possible stable 
particles, whatever the primary process. For instance, even if a tree level annihilation into electron/positron is considered, in 
the final spectrum also antiprotons will be present. Quantitatively, at the TeV scale and beyond EW corrections of infrared 
origin are typically as big as the tree level values and cannot therefore be neglected. These corrections have been considered in 
the context of DM signals [121 1131 [14] : however corrections growing like s/m%y were found, while we find corrections featuring 
double logarithmic growth. The reasons for these discrepancies are analyzed in Sec. [4] 

The purpose of this paper is to investigate the impact of EW radiative corrections for possible TeV scale DM signals, namely, 
trying to contribute to give an answer to the question raised above. Since we focus on the impact of EW corrections, we do 
not aim at constructing a realistic DM model, nor do we consider the effects of propagation from the interaction point to the 
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detection point. Rather, we consider a very simple model: a heavy Z' gauge boson corresponding to a U'(l) group factorized 
with respect to the SM group. We only consider Z' decay into a fermion-antifermion pair, to which we add a weak gauge boson 
emission and one loop radiative corrections. 



2 Gauge boson emission in the soft /collinear region 

We add to the standard model Lagrangian Csm a vector boson Z' with mass M bigger than 1 TeV, belonging to an extra U'(l) 
gauge symmetry and singlet under the SU(3)c ® SU(2)_l C8> U(l)y gauge symmetry. The relevant couplings with quarks and 
leptons are dictated by the property of gauge invariance. Given the usual SU(2)_l doublet L = (vl, &l) T and the singlets en, vr 
, we have 

Cint. = Z'^J£- J£ = fiZ^L + f&SrfvR + f&R^en, (1) 

with similar expressions holding for the other families and for quarks. Let us consider the case fh = fit = Oj so that the Z' 
couples to left electron and neutrino with equal strength Jl ■ 

We indicate with T2 the width for the process Z' — > ee; we wish to calculate the effect of adding one weak gauge boson 
emission. As is well known, in the high energy regime M S> mw the leading contributions to the three-body width Z'(P) — > 
e + (pi)W~ {k)v(p2) are produced by the region of the phase space where the emitted boson is collinear either to the final fermion 
or to the final antifermion; moreover the three-body width is factorized with respect to the two-body one in this region. Here 
we show explicitly this factorization and we calculate the expression for the three-body width. 

Let us show that the phase space factorizes in the region where the gauge boson momentum k is collinear to the emitting 
fermion momentum p\ (the region where k is collinear to P2 can be treated in the same way). In fact in this region p\ — k 2 = 
implies p 2 — (masses neglected); moreover we can write 

Mpj(M^ 5(4){p - (£h 5(i)[p - p - p2) w^L (ih) - v 

where P is the decaying Z' momentum and x is the fraction of energy carried away by the gauge boson, ko = xPo- The 
three-body phase space therefore factorizes with respect to the two-body one in the collinear regio: 
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'(27r) 3 2fc° Vl-a; 

Furthermore, the amplitude squared factorizes as well, as we will show now 




Figure 1: Amplitude squared that contributes to the primary particle spectrum (see text), (a)-(c) describe real 
emission contributions, while (d) is an example of a virtual correction. 



Let us first consider the contribution to the modulus squared of the amplitude shown in Fig. [Tp. This contribution can be 
written as 

where D va — —g va + k v k t7 /m 2 v is the sum over the emitted W physical polarizations and e (J (P) is the physical Z' polarization. 
In the following we will systematically neglect terms that, when integrated over the phase space, produce contributions not 

iFor convenience, we include a 1/(2M) factor in the definition of dT: dT n = l/(2M)\M n \ 2 d& n , with d<J> n being the usual phase space. 
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growing with energy: the symbol ~ refers to this approximation. Let us first consider the contribution coming from the g va 
component of the sum over polarizations. In the collinear approximation k ~ x/(l — x) pi after some Dirac algebra we obtain 

- M^^) ^- p2Tr < (5) 

which features factorization in the collinear region, since because of (J3]): 

^fc-.,- - (itt [{k+PlY - + Pa)a - my ih {/ V^WWTr [^ h YP L ] } , (6) 

and the integral between braces in Eq. © is precisely the tree level width r 2 . In the collinear/infrared region we have (see the 
Appendix) : 

d 3 Je 2pi -p a dx 1 , M 2 x 2 

* " " 2 ln 7^> ( 7 ) 



(27r) 3 2w [(pi + k) 2 - raj] [(p 2 + fc) 2 - ra 2 ] x 16tt 2 4m 2 , 

so that we finally obtain 



dF& __„ ^F 2 ——U(l-x)]nZ-^-\. (8) 
where aw = 5 1 /4ir. Analogous calculations produce the contribution coming from the k^kv term: 

The remaining contributions represented in Figs. [TJa) and[TJb) are given by 

dr f _ 9 ^r 2 ^^( 2a; 2 in^-l, (io) 

,„a i — i aw da; f„ , , M 2 x 2 M 2 x 2 \ 

dFf _ k k ,2 «T 2 -f- I 2a;(l-a;)ln— + -r— 5- \, 11 

and 

, r s _ "w f M 2 x 2 ) 

dr E poi =w™? v ~ F2 ^tt \ 4<r; • (13) 

Finally, the sum of all amplitudes squared, generated by Figs. [Ua)-[ljc), gives 



aw dx ( . , x 2 M 2 



Notice that, although terms proportional to M 2 /rriw are present in the single contributions ((81 ll3p . these terms disappear from 
the final result (fli)) . 



3 Virtual corrections and primary particle spectrum 

The calculation of virtual corrections, which must be included in order to predict the primary particle spectrum, poses no 
particular difficulty. In the very high energy regime we are considering, they are dominated by the region of integration over 
the virtual momentum where the exchanged gauge boson [see Fig. QJd)] is close to the mass shell and has the same kinematical 
structure of real emission. Virtual corrections are thus dominated by the soft/collinear region just as the real emission ones 
and factorized with respect to the tree level amplitude. Referring the reader to the relevant literature [151 116j . we wish to point 
out that, because of (a) factorization and (b) unitarity of the theory, virtual contributions can be derived from real emission 
calculations described in the previous section. 

The fully inclusive decay width Ttot(Z' — > ff + X) is given, at the order of perturbation theory considered here, by the sum 
of the width for the case of no emission \To(Z' — > //)] and the one with one gauge boson emitted [Ti(Z' — > ff+X, X = 7, Z, W)\. 
Because of factorization in the leading collinear/infrared regime, one can write Yi = YBomPi, where TBom is the tree level 
value and Pi are functions of couplings and energy scales. Now, since the theory is unitary, one has 

TtOT = Ffl orn (Po + Pi) = rgorn, (15) 
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so that the inclusive cross section equals the tree level value and Pi can be interpreted as probabilities. Now, Pi can be found 
by integrating Pr(x) in Eq. (|14[) over the available phase space and 



Pi = 1 



dxP™ (x) : 



1 



aw 
4% 



1 2 M 

2 Ami, 



3 M 
2 ln S 



(16) 



If we indicate with z = 1 — x the momentum fraction carried away by the positron, virtual corrections are described by a 
distribution peaked at z = 1, so that one finally obtains 



aw 
4-k 



M 



1 l 2 

2 ln W" 



3 M 

2 4^: 



(17) 



J^(«) = fl,«(l-«) = *(1-«)^1- 

' w * 

The spectrum of the positron is described by the distribution P w (z) = Pr(z) + Pv{z). To this one must add the effects coming 
from radiation of a Z boson, which can be derived in a similar way, and the effects of photon radiation, that we discuss below. 

In the case of QED the kinematic is different, since the photon is massless and since collinear singularities are cut off by the 
emitting particle mass (m/, considered here to be a fermion). The distribution of emitted photons, derived in the Appendix, is 



aW 2tt Am) 
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(18) 



where a = e 2 /Atv and where e is an infrared regulator, with dimensions of a mass, having the physical meaning of lowest energy 
for the photon. The distribution for virtual corrections can be derived using (|15[) and is 
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The QED distributions depend on the arbitrary parameter e and are divergent in the limit e — > 0, so they are obviously 
unphysical. As is well known, the way out is to introduce a finite resolution AE on the observed hard particle (say, a positron). 
The physical meaning is that what is really observed is not a positron alone, but rather a positron together with a soft photon 
of energy u < AE. We take this into account by substituting P v + P^ in the region z > 1 — "^jf- with a flat distribution whose 
integral is the same as the one of P v (z) + PZ{z) in that region: 



a M 2 

( z ) = 7T- m ~, — ? 

2tt Amj 
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Let us now discuss which value one should give to AE. In a collider, this value would be given by the (known) characteristics 
of the detector. Here however, the situation is less clear because of the effects of propagation from the interaction region to the 
point where the detector is physically placed. We choose to analyze two values for AE. The first could be called the "optimal 
resolution case": AE is chosen in such a way that the distribution (|20[) becomes a continuous line; in the case at hand this 
corresponds to AE w 0.3 MeV0. The resulting distribution for the positron is drawn in Fig. [2] where contributions from 7, 
W and Z radiation have been added together. Since the actual resolution on the positron energy is certainly worse than the 
"optimal" one, it is possible to obtain the actual distribution from the one in Fig. [2] simply by dividing it into "bins" of finite 
width. This is done in Fig. [3] where the more realistic case AE— 30 GeV has been chosen. Finally, the region z ~ 1 from 
the "ideal case" in Fig. [2] is drawn for convenience in Fig. [4] here the contributions from 7, Z and W radiation are drawn 
separately. 

Let us now discuss our results. In the first place, it is apparent from Fig. [5] that the emission of weak gauge bosons plays 
a significant role in determining the spectrum of the primary particle; therefore one should always consider QED (7) and weak 
(W,Z) radiation together at very high energies. This is true even more so if the primary particle is an EW gauge boson itself 
instead of a fermion. In fact in this case QED radiation is partially suppressed because the gauge boson mass provides the 
collinear cutoff in Eq. (|18[) . in place of the much smaller fermion mass considered here. Moreover, in the case of final EW gauge 
bosons since EW corrections of infrared origin are proportional to the Casimir of the external legs representations [S], they are 
expected to be bigger in magnitude. 

We can see from Fig. [2] that even after inclusion of EW radiative corrections, the antifermion spectrum is rather sharply 
peaked at z = 1. This can be seen also from Fig. [3] where the tree level distribution (dashed line) is also plotted for convenience. 
Virtual EW corrections deplete the first bin (from the right) by about 30 %, and produce therefore a significant effect. However 
the second bin is depressed with respect to the first one by 1-order of magnitude, and the others are even lower, so that the 
great majority of events falls into the first bin even after including radiative corrections. 



Figs. I2l4l3l are plotted in the case of the spectrum of an antimuon and a DM mass M = 10 TeV. 
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4 Comparison with existing calculations 



In |12] the effect of adding a weak gauge boson emission to a DM annihilation cross section at very high energies y/s ^> mw 
was considered; in [14] a similar calculation for the case of very heavy decaying DM [M 3> mw) wa s done. In both cases, 
corrections growing like the square of the cm. energy were obtained. In particular, it was found that the cross section (width) 
with one additional gauge boson in the final state is obtained by multiplying the original one by a factor Rz,w(M): 

M 2 

Rz,w(M) = a w K z ,w —2- + -, (21) 
m w 

where M is the relevant high energy scale (y/s in the case of scattering, DM mass in the case of decay), is supposed to be much 
higher than the weak scale; Kz.w is a constant that depends on whether a Z or a W is radiated. The dots here stand for terms 
that are subleading in the M 2> mw regime. 

On the other hand, when integrated over the final gauge boson phase space (variable a;), Eq. (|14[1 gives 



Rw(M) - — 



M 2 , M 

3 In —^5- + . 



4771^ 4m 2^ 



(22) 



r 2 (Z' -> ee) v ' 8tt 

Similarly, virtual corrections described by (|17|) grow like the square of the logarithm of M 2 /m w . Our result therefore disagrees 
with the results obtained in |12l 1131 fl4] . where corrections growing like M 2 /m w were found. Here we wish to point out that 
the discrepancy is related to the introduction of dimension-4 operators that break SU(2)l <g> U(l)y gauge invariance, and that 
one runs into ambiguous and possibly inconsistent results when trying to calculate virtual corrections in such framework. 

Terms growing like M 2 /m^ are indeed present in the single contributions to the amplitude squared [see (|9llip ] and are 
related to the terms in k^k^/m^ in the sum over the emitted gauge boson polarizations. However, such terms are absent from 
the final result (|22[) . This is a consequence of gauge symmetry in the form of Ward identities that are depicted in Fig. [5] These 
identities relate on shell amplitudes with an external gauge boson with corresponding amplitudes with an external Goldstone, 
in the following way: 

■^M l *(k,...)=iM( V > w (k),...), (23) 
mw 

with mw being the mass of the relevant (W or Z) gauge boson. 

Since Goldstone bosons couple with fermions through their mass, the right-hand side is close to zero. Therefore the terms in 
the polarization sum proportional to fc M fc„, that are formally dominant by power counting with respect to the g^ v term, are 
strongly suppressed at high energies. 

A common feature of [121 114] is the introduction in the Lagrangian of dimension-4 operators that explicitly break SU(2)l ® 
U(l)y gauge invariance: D^Vl^^vl and/or D(Vrvl + VEvr) [to be compared with our gauge symmetry invariant Lagrangian 
in 0]. Let us first consider real emissions. 

Suppose that, following [TJ], we introduce a D^VZ^fz, interaction in the Lagrangian without its electron counterpart. 
Clearly, Ward identities are broken: for instance, from the point of view of our gauge invariant example, in the diagrams 
of Fig. [5] the second one on the left-hand side is missing. Then, the terms proportional to M 2 /m w generated by the term 
proportional to fc M fc„ in the sum over polarizations do not cancel between the various contributions and are eventually present in 
the expression for the width. What one is really doing here is largely overestimating the contribution from longitudinal degrees 
of freedom, whose polarization is = + O (^jjjp) and whose contribution, although leading by naive power counting with 
respect to transverse, is suppressed in a gauge invariant theory because of Ward identities. 

As we have seen , virtual corrections have to be considered together with real emission for calculating the observed spectrum 
and, by virtue of (|15[) . are important for unitarity of the theory. However we wish to point out that if, as done in [13], one 
calculates virtual corrections, the consistency of the calculation itself is put into doubt because, since Ward identities are broken, 
the final results will depend on the (arbitrary) gauge that one chooses in order to give a prescription for the weak gauge boson 
propagators. In fact the infrared structure is dominated by the region of phase space with on-shell gauge bosons. Since Ward 
identities are broken, the gauge dependent terms oc fc M fc„ in the propagators are completely different, say, in the Feynman gauge 
and in the unitary gauge. These problems are generated by the exchange of soft quanta with energies greater than, but close 
to, the weak scale: therefore, we think that they cannot by cured by any "UV completion" beyond the TeV scale. 

In Ref. P3], for phenomenological purposes, a chirality- (and gauge symmetry-Wiolating operator of the form D(J7rvl + 
Vlvr), with D being a scalar and gauge singlet, is introduced into the LagrangiarU- In this case only tree level gauge boson 
emissions have been considered. Notice that it is true that radiative corrections generate chirality violating terms proportional 
to the gauge symmetry breaking vacuum expectation value v through fermion masses [JJ However, since this term is generated 
in the framework of a gauge invariant theory (like the SM itself), Ward identities are obviously respected. On the contrary, 
by introducing a dimension-4 operator that explicitly violates gauge invariance, Ward identities are broken and by trying to 
compute radiative corrections one runs into the problems signaled above. One can, of course, write a gauge invariant interaction 



^In a gauge invariant context 1171 . a milder growth with energy of contributions coming from the emission of a gauge boson was found. 
"For some surprising features related to the high energy behavior of such terms we refer the reader to |18| . 
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if the scalar DM particle is an isospin doublet: this would be similar to the usual Higgs-fermion-antifermion coupling in the 
standard model. We wish to point out that in this case, like in any other case where the DM particle has a weak charge, gauge 
boson emissions from the initial state have to be considered together with emissions from the final legs (see Fig. [6}. If only final 
state radiation is considered, Ward identities are broken and one effectively breaks the SU(2)^ <g> U(l)y symmetry, potentially 
generating again "spurious" terms growing like the square of the energy. 

One final comment: in our case the Z' can perfectly well couple to right neutrinos only, with the choice fz. = /J? = 0, 7^ 0. 
The Z' then decays only into neutrinos and this situation, barring possible tiny chirality breaking effects proportional to fermion 
masses, is left unchanged by radiative corrections. Therefore, no final electrically charged states are present, even after including 
radiative corrections. 

5 Conclusions 

In this work we have analyzed the impact of radiative electroweak corrections on the spectrum of the final products resulting 
from the decay of a very heavy (M > I TeV) weakly interacting particle. Determining accurately this spectrum is an important 
issue, in view of recent experimental results that can be interpreted as a dark matter signal. We have considered a simple 
model with a Z' gauge boson decaying into leptons, and rediscussed one loop radiative corrections plus emission of a weak 
gauge boson. We have found that electroweak corrections play a relevant role in this game; more precisely we can summarize 
our conclusions as follows: 

• EW radiative corrections at one loop are of the order of 30 % in the considered case of fermions as primary particles, 
and grow like the log squared of the DM mass. One expects these corrections to be even bigger in the case of EW gauge 
bosons as primary particles (see Sec. |3); in any case, higher order effects play a significant role. 

• QED corrections and "pure weak" corrections produced by W and Z exchange have a similar impact on the spectrum of 
the primary particle, so they should always be considered together. 

• Different from recent results in the literature, we do not find corrections growing like the square of the energy. We have 
shown that the latter result can only be obtained if the Ward identities related to SU(2)l CS> U(l)y gauge symmetry are 
broken. We have argued that this choice can lead to ambiguous results for virtual corrections and we think that this 
point should be better clarified in the future in view of its importance. 

• As already noticed in |14] . since EW corrections link all SM particles, all stable particles (including antiprotons) will be 
present in the final spectrum, independently of the primary particles (say, lepton/antileptons) initially considered. 

• Heavy DM annihilation/decay can produce final states in which no electrically charged particles are present (see Sec. 2}. 

Let us now comment on the possible impact of EW corrections on indirect signals of dark matter. A detailed quantitative 
analysis including the effects of propagation goes beyond the scope of this work, but we can make some general comments. 
The first point is related to the fact that the e excess measured in recent data in the 10 - 100 GeV region [6j [9] is not 
accompanied by a corresponding excess in the antiproton flux in the same energy region: the fraction 2 is indeed compatible 
with the astrophysical backgrounds [5J. One could then consider for instance DM annihilating to leptons only, or to W + W~ 
pairs with so high an invariant mass that the antiproton excess is shifted to energies higher than 100 GeV [19!. In both cases, 
however, taking into account EW corrections changes the scenario: in the first case allowing for final antiprotons coming from 
the decay of the weak gauge bosons irradiated by the leptons, and in the second case by softening the final gauge bosons 
spectrum. So radiative corrections of EW origin should be definitely taken in such kind of analyses. 

Another issue is whether it is necessary, for the problem at hand, to consider higher order EW corrections. At first sight, 
one loop virtual corrections plus real gauge boson emission, as are done here, seem to be sufficient. In fact, as we have seen, 
one loop corrections are of the order of 30 % and higher order virtual corrections should not change the picture much: the 
experimental precision of the measurements involved is much lower than that of a earth based accelerator of course. However, 
this conclusion turns out to be a hasty one. In fact we are interested in the particles' spectra, in the experimentally interesting 
region 10-100 GeV, produced by a decay annihilation taking place at a high invariant mass of the order of 1 TeV or more. 
Then, even if EW corrections on, say, the total cross section are relatively small, the impact on the spectrum at low energies 
can be dramatic, since "soft" (~ 10 GeV) particles are radiated copiously in the presence of a "hard "(~ TeV) process. So one 
can determine whether or not higher order EW corrections are relevant in this context only after calculating their effects: this 
can and should be done in the future. Higher order EW radiative corrections and resummation of leading effects in the high 
energy regime have been extensively studied in recent years [31 E]i recently a way of systematizing virtual corrections through 
soft collinear effective theory has been considered [20] . However, EW corrections are usually studied in two rather "extreme" 
cases: virtual corrections to a given hard process (in this case, Z' decay into a hard fermion and a hard antifermion) or the 
inclusive case where soft weak gauge boson emission is fully included. Determining the spectrum resulting from heavy DM 
annihilation/decay is a different, and difficult, case, where one wishes to know the distribution of a large number of particles 
in the final state. One possibility to deal with such a problem is to implement the EW evolution equations [IB] (the analogous 
of the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi equations in QCD) in parton shower Monte Carlo models. This seems to 
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be unavoidable in the nearby future, especially if high energy lepton colliders [2T], where showers generated by electroweakly 
charged particles will be ubiquitous, come closer to see the light. 



6 Appendix: Kinematics 

We consider the decay Z'(P) -t /(pi)/(p2) with P = (M, 0), pi = (M/2, p), p 2 = (M/2, -p), and p 2 >2 = mj. Then we consider 
the additional emission of a gauge boson with momentum fc and mass A. The latter is parametrized in the following way: 

k = (cj,k), fc 2 = A 2 ; p ■ k = \p\ fc cos 6; (24) 

We consider the region where the emitted boson is collinear to p 1 = p so that 9< 1 (the region where A; is collinear to p 2 = —p 
can be treated in a similar way). In our calculations the following expression appears: 

d A ~£ ggi ■ p 2 , . 

aT > ( zo ) 



(2tt)32l; [(pi + k) 2 - m}][(p 2 + k) 2 - mj] 
where the differential phase space is evaluated in the collinear region to give 



d 3 k u 



(2tt) 3 2u; (4tt) ; 



■dude 2 . (26) 



The denominator appearing in (|25p is evaluated differently in the case of QED (A = 0) and in the case of weak gauge bosons 
emission (A = mw 3> rrif). 

In the case of QED we have 

(pa + k) 2 - m 2 f = 2p 2 ■ k w 2Muj, (p a + fc) 2 - m) = 2 Pl ■ fc w ^ (0 2 + ^), (27) 



and the integration over # 2 readily produces 



d 3 k 2pi ■ p 2 = dx _J_ M 2 

(2tt) 3 2o; [(pi + k) 2 - m 2 f ][(p 2 + k) 2 - mj] ~ x 16n 2 * Am 2 ' ( ' 

where x = 2ui/M is the fraction of energy carried away by the emitted photon. 
In the case of weak gauge boson emission we have 

(pa + k) 2 - m) = m 2 w + 2p 2 • k w 2Mu, (pi + fc) 2 - m) = 2pi ■ fc + w ^(6> 2 + ^), (29) 



and integrating over : 



d 3 k 2pi ■ p 2 = m2x2 (o^ 

(2ti)^2uj [(pi + fc) 2 - m 2 ][(p 2 + k) 2 ~ m 2 } x 16tt 2 n Am 2 w ' 1 j 
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Figure 2: Spectrum ^P- of the primary particle (antimuon) coming from the decay of a heavy Z' [M — 10 TeV) after 
inclusion of one loop virtual corrections and EW gauge bosons emission. Here and in the following figures i— 2E/M, 
with E being the antimuon energy. An "optimal resolution" is assumed (see text). 
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Figure 3: Continuous line: the same as Fig. [2] but with an experimental resolution AE(/i) — 30 GeV. Dashed line: 
tree level distribution with no electroweak corrections. 
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Figure 4: The same as Fig. [5] but zoomed in the region z ss 1 and with the contributions coming from 7, W and Z 
emission drawn separately. 
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Figure 5: Ward identity relevant for the process considered in this paper. 
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Figure 6: Gauge boson emission diagrams (the wavy line is a Z or a W) for the case of a heavy scalar decaying into 
a fermion or an antifermion. 
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